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Abstract 

Background: Microglia, the resident immune cells of the central nervous system (CNS), have two distinct 
phenotypes in the developing brain: amoeboid form, known to be amoeboid microglial cells (AMC) and ramified 
form, known to be ramified microglial cells (RMC). The AMC are characterized by being proliferative, phagocytic and 
migratory whereas the RMC are quiescent and exhibit a slow turnover rate. The AMC transform into RMC with 
advancing age, and this transformation is indicative of the gradual shift in the microglial functions. Both AMC and 
RMC respond to CNS inflammation, and they become hypertrophic when activated by trauma, infection or 
neurodegenerative stimuli. The molecular mechanisms and functional significance of morphological transformation 
of microglia during normal development and in disease conditions is not clear. It is hypothesized that AMC and 
RMC are functionally regulated by a specific set of genes encoding various signaling molecules and transcription 
factors. 

Results: To address this, we carried out cDNA microarray analysis using lectin-labeled AMC and RMC isolated from 
frozen tissue sections of the corpus callosum of 5-day and 4-week old rat brain respectively, by laser capture 
microdissection. The global gene expression profiles of both microglial phenotypes were compared and the 
differentially expressed genes in AMC and RMC were clustered based on their functional annotations. This genome 
wide comparative analysis identified genes that are specific to AMC and RMC. 

Conclusions: The novel and specific molecules identified from the trancriptome explains the quiescent state 
functioning of microglia in its two distinct morphological states. 



Background 

Microglia are the prime immune effector cells of the 
central nervous system (CNS). The origin, morphology 
and role of microglia in health and disease were first 
elaborately described in 1939 [1]. The amoeboid micro- 
glial cells (AMC), which are abundant in the periven- 
tricular white matter, namely the corpus callosum (CC) 
of the brain function as macrophages in the developing 
brain. Studies have demonstrated that AMC gradually 
transform into ramified microglial cells (RMC) with ad- 
vancing age [2]. 
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The time course of development of microglia differs in 
different regions of the brain [3]. In the CC, AMC pre- 
ponderate a week before birth in mice and rats [4-6] and 
actively phagocytose the cellular debris and refine axonal 
connectivity during the first postnatal week [7-9]. This is 
followed by the gradual transition into RMC, which sur- 
vey the brain parenchyma with their fine and non- 
overlapping ramifications, thereby monitoring chronic 
and acute insults [10]. Upon activation by trauma, infec- 
tion or any other neurodegenerative stimuli, microglia 
retract their ramifications, transform into amoeboid or 
spherical shape, produce pro -inflammatory cytokines 
and display phagocytosis [11]. These microglia are 
known to be activated or reactive. 

Studies in chronic, aging associated neuropathologies 
such as Alzheimer's disease (AD) [12-14], and Parkinsons 
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disease (PD) [15] indicate persistent microglial activation as 
the major causative factor in disease exacerbation. Aging 
brains are often characterized by the presence of primed 
microglia, which present an altered cytokine profile in com- 
parison to their counterparts in younger brains [16,17]. 
These microglia produce an exaggerated inflammatory re- 
sponse when activated [18,19] leading to prolonged cycles 
of proliferation and production of pro-inflammatory 
cytokines which eventually render them neurotoxic. Fur- 
ther, chronic microglial activation has been shown to cause 
the impairment of adult neurogenesis in hippocampus [20] 
and damage to the periventricular white matter (PWMD) 
in the early postnatal brain [21]. Hence activated microglia 
in both postnatal and adult stages can have neurotoxic 
effects on the CNS by causing excessive inflammation. 
Identification of ways to attenuate microglia-mediated 
neuroinflammation, therefore, has been the primary consid- 
eration in therapeutic strategy. There is accumulated infor- 
mation on the factors that contribute to the activation, 
migration, proliferation and immune response of microglia 
over the years [22,23], but the gene expression and signal- 
ing networks that function within these cells are yet to be 
fully clarified. 

Gene expression profiles of microglia from primary cul- 
tures are available, but their expression profiles have been 
found to be altered once isolated from their natural milieu 
[24]. It is striking that investigation on the expression pro- 
files of functioning genes of AMC and RMC in vivo in their 
quiescent state have remained elusive. In this connection, 
we carried out a global gene expression profiling of AMC 
and RMC in situ by isolating them from the CC of rat brain 
using laser-capture microdissection (LCM) tool. Overlap- 
ping the transcriptome onto several online and commercial 
databases, our current study aimed to identify molecular 
candidates that are associated with the morphological 
transformation and physiological functioning of microglia 
in the developing brain. Further, we have identified the 
genes that render sternness' and monocytic functions to 
AMC and RMC. The transcriptome profiling has also led 
to the identification of several genes that may be vital in 
regulating microglial proliferation, differentiation, migra- 
tion, and ramification. 

Methods 

Ethics statement 

In the handling and care of animals, the International Guid- 
ing Principles for Animals Research, as adopted by the 
Institutional Animal Care and Use Committee (IACUC), 
National University of Singapore, were followed. All efforts 
were made to minimize pain and the number of rats used. 

Laser-capture microdissection (LCM) 

Whole fresh brains were removed from 5 -day postnatal 
Wistar rat pups (n = 3) and 4-week old Wistar rats 



(n = 3) and placed in liquid nitrogen immediately for a 
short time and then frozen in a cryostat (Model No. CM 
3050 S, Leica Microsystems GmbH, Wetzlar, Germany). 
The forebrain was sectioned coronally through the CC 
at 5 (im thickness and mounted on precleaned slides. 
The sections were fixed in 75% ethanol for 1 min and 
incubated with peroxidase conjugated isolectin (1: 50, 
Cat. No. L5391, Sigma-Aldrich Co., MO,USA) for 
15 min. The sections were then dehydrated by a graded 
series of ethanol and cleaned in xylene. The slide was 
placed on the microscope stage of MMI CellCut (Mo- 
lecular Machines & Industries, Glattbrugg, Switzerland). 
The 4 X, 10 X to 40 X objective lenses were used to 
achieve the proper placement of the cap (for cell collec- 
tion) above the CC. Lectin stained microglia cells (AMC 
from 5 -day and RMC from 4- week old rat brain CC) 
were selected and cut by laser and collected into the cap 
of tube (Cat No. 50202, Molecular Machines &Indus- 
tries, Glattbrugg, Switzerland). Extra care was taken to 
minimize the contamination of materials from other cell 
types while laser dissecting microglia from the CC. 

Microarray analysis 

Total RNA was extracted from 600 isolated microglia 
cells per group using RNeasy micro kit (Cat. No. 74004, 
Qiagen, CA, USA), quantified by Nanodrop 1000 
(Thermo Scientific, MA, USA) and hybridized to each 
microarray chip. RNA (15 ng) was reverse transcripted 
into the first-strand cDNA using a T7-OHgo (dT) Primer 
(Two-Cycle Target Labeling and Control Reagent pack- 
age, Affymetrix, CA, USA). After second-strand cDNA 
synthesis, the double-stranded cDNA was purified and 
served as a template in the first cycle of in vitro tran- 
scription (IVT) reaction. The unlabeled cRNA was then 
reverse transcripted into the first-strand cDNA of the 
second cycle using random primers. Subsequently, the 
T7-OHgo(dT) Promoter Primer was used in the second- 
strand cDNA synthesis to generate double-stranded 
cDNA template containing T7 promoter sequences. 
Then the double- stranded cDNA was amplified and 
labeled using a biotinylated nucleotide analog/ribonu- 
cleotide mix in the second IVT reaction. The labeled 
cRNA was then cleaned up, fragmented, and hybridized 
to Rat Genome 230 2.0 Array (Cat. No.900506, Affyme- 
trix, CA, USA). A total of six arrays (three each for 
AMC and RMC) were carried out in the present study. 
The arrays were stained according to the manufacturer s 
protocols and then scanned with the Genechip scanner 
(Affymetrix, CA, USA). Initial analysis of the scanned 
images was performed by GeneChip Operating Software 
(GCOS, Affymetrix, CA, USA). For absolute analysis, 
each chip was normalized to a target intensity of 500, 
and probe sets were assigned a signal intensity and de- 
tection call of Present, Marginal or Absent. 
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Data analysis and generation of gene lists 

The absolute data (signal intensity, detection call and de- 
tection p-value) were exported into GeneSpring GX 7.3 
software (Agilent Technologies, CA, USA). All the six 
chips were globally normalized and the genes of over 
2-fold differential expression were filtered out and used 
for functional analysis. 

Data normalization and generation of gene lists using 
MATLAB 

Raw CEL files of the six chips were RMA (Robust Multi- 
chip Average) normalized using the Affymetrix Expres- 
sion Console Version 1.1 (Affymetrix, CA, USA). The 
normalized data was then used to identify differentially 
expressed genes between AMC and RMC in MATLAB 
R2009a (MathWorks, MA, USA). For the statistical ana- 
lysis, we used the 'Exploring Gene Expression Data 
demo scripts in the Bioinformatics Toolbox . The data 
was filtered for removing genes with low expression 
values and low variance across chips. Further, t-test was 
performed to retain genes with p-values less than 0.05 
and a Volcano Plot was generated to identify the two- 
fold differentially expressed genes. The microarray data 
discussed in this publication is MIAME compliant and 
has been deposited in NCBIs Gene Expression Omnibus 
(GEO, http://www.ncbi.nlm.nih.gov/geo/). It is accessible 
through GEO Series accession number GSE29885. 

Gene expression profile clustering and pathway analysis 

Agglomerative average-linkage hierarchical clustering of the 
different experimental groups was obtained for selected 
groups of genes with GeneSpring GX 7.3 software (Agilent 
Technologies, CA, USA) with standard correlation used as 
the similarity matrix. The gene lists obtained was fed into 
Pathway Studio 6 software (Ariadne, MD, USA) to generate 
pathways for identifying interactions between the genes for 
validation purposes. 

Analysis of gene lists 

The gene list generated from MATLAB was used to 
identify functional groups enriched in the AMC and 
RMC using DAVID Bioinformatics Database [25,26]. To 
identify the 'Sternness' of AMC and RMC, we compared 
our gene lists to gene lists enriched in embryonic, neural 
and hematopoietic stem cells [27]. Since the data were 
accumulated from a different microarray platform, we 
found orthologs to their genes pertaining to our plat- 
form using the online NetAffyx application (Affymetrix, 
CA, USA). For comparison of our gene expression data 
to that of peripheral blood monocytes [28], the raw CEL 
files of monocyte expression data were downloaded from 
NCBI GEO (Gene Expression Omnibus) and the ortho- 
logs pertaining to our platform were identified using the 
online NetAffyx application. These files were RMA 



normalized in Affymetrix Expression Console Version 
1.1 (Affymetrix, CA, USA) and subsequently the average 
expression values of the monocyte genes were compared 
to our microglia gene lists. 

Double immunofluorescence staining on postnatal rat 
brain sections 

5 -day and 4- week old Wistar rat pups were purchased 
from the Laboratory Animal Centre, National University 
of Singapore. The animals were perfused and fixed with 
4% paraformaldehyde for further procedure. For double 
immunofluorescence staining, forebrain sections at 30um 
were cut through the corpus callosum using cryostat 
(Model No. CM 3050 S, Leica Microsystems GmbH, 
Wetzlar, Germany). The sections were incubated with 
purified mouse anti-OX-42 Ig (1:50; Cat No. CBL1512, 
Millipore, MA,USA) along with rabbit anti-ETO (1:100; 
Cat No. sc-28693, Santa Cruz Biotechnology, Inc. CA, 
USA) or with rabbit anti-Dcx (1:100, Cat No. abl8723, 
abeam, Cambridge, UK) or with rabbit anti-Sox4 (1:100; 
Cat No.sc-20090,Santa Cruz Biotechnology, Inc. CA, 
USA) or with rabbit anti-Soxll (1:100; Cat No.sc-20096 , 
Santa Cruz Biotechnology, Inc. CA, USA) or with rabbit 
anti-Sept9 (1:100; Cat No. sc-130263, Santa Cruz Biotech- 
nology, Inc. CA, USA) with rabbit anti-Sept4 (1:100; Cat 
No. sc-20179, Santa Cruz Biotechnology, Inc. CA, USA) 
overnight at 4°C. On the following day, the sections were 
further incubated with either FITC-conjugated goat-anti- 
mouse IgG (1:100; Cat No. F9137, Sigma- Aldrich Co., MO, 
USA) or Cy3-conjugated sheep-anti-rabbit IgG secondary 
antibody (1:100; Cat No. C2306, Sigma-Aldrich Co, MO, 
USA). The sections were counterstained with DAPI (1 ug/ 
ml, Cat. No. D1306, Invitrogen, CA, USA) and mounted 
with a fluorescent mounting medium (DakoCytomation, 
Glostrop, Denmark). Photo-images were captured using a 
confocal microscope (Olympus FV1000, Tokyo, Japan). 

Cell culture 

BV-2 cells (a widely used murine microglial cell line) 
were maintained at 75 cm 2 culture flasks in Dulbeccos 
Modified Eagles Medium (DMEM, Sigma, St. Louis, 
MO, USA; Cat. No. 1152) supplemented with 10% fetal 
bovine serum (FBS, HyClone, Logan, UT) and cultured 
in 37°C in a humidified atmosphere of 5% C0 2 and 95% 
air incubator. Cells were seeded on 6-well plates at about 
a density of 1.0 x 10 6 per well for RNA isolation. 

RNA isolation and real-time RT-PCR for validation of 
microarray data 

Total RNA from laser-captured AMC and RMC was 
extracted using miRNeasy Mini Kit (Qiagen, Germany, Cat. 
No.217004) and RNA from BV-2 cells was extracted with 
RNeasy Mini Kit (Qiagen, Germany, Cat. No. 75161) 
according to the manufacturer s instructions and quantified 
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Figure 1 A-F. Identification and isolation of amoeboid microglial cells (AMC) and ramified microglial cells (RMC) from the corpus 
callosum (CC) of 5-day old and 4-weeks old rat brain respectively. Figure A shows AMC and Figure D shows RMC stained with lectin under 
laser capture microscopy. These cells are laser-cut along their periphery and isolated. Figure B and E show the region of the stained tissue section 
wherein the cells have been removed, and figure C and F show the isolated cells collected in the cap of vial. Arrows indicate the same cells in all 
the three images. Scale bars: A-F 50 urn. G. Correlation plots. Correlation plots were generated in Affymetrix Expression Console 1.1 after RMA 
normalizing raw CEL files of AMC and RMC expression data. The color scale indicates the degree of correlation between two different samples. A 
value of close to 1 refers to a high correlation. 



spectrophotometrically. 2 \ig of RNA from each sample was 
added to a total volume of 25 ul reaction mixture contain- 
ing 2.5 uM of oligo (dT) primer (Promega, Madison, WI 
USA; Cat No. C110A), and 200U of Molony Murine 
Leukemia Virus Reverse Transcriptase (M-MLV, Promega, 
Madison, WI, USA; Cat. No. M5314). The reaction was 
initiated by incubating the reaction mixture for 1 h at 42°C 
for reverse transcription, and stopped by heating for 
10 min at 70°C. Aliquot (0.5 ul) of the each reverse tran- 
scription product was added to the 10 ul reaction mixture 



containing QuantiTect SYBR Green I, 0.5 uM of each pri- 
mer corresponding to Runxltl, Sept9, Sept4, Mbp (rat), 
Gapdh, Dcx, Mbp (mouse), or (3-actin and 4 mM MgCl 2 to 
amplify the genes in ABi 7900HT Fast PCR system (Ap- 
plied Biosystems, USA). The primer sequences of Runxltl 
are forward: 5 -ACGAACAGCTGCTTCTGGAT-3 and 
reverse: 5 -TGCTTGGATGTTCTGAGTGC-3 ', Sept 9 are 
forward: 5 -AACCATGTCCCTCGAACTTG-3' and re- 
verse: 5 -AAGAGAGAGGGGACACGACA-3 ', Sept 4 
are forward: 5 -CTCATCCGGGAGAAAGATGA-3 ' and 
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Figure 2 A. Cluster Analysis. Cluster analysis shows changes in gene expression profiles of AMC and RMC Agglomerative average-linkage 

hierarchical clustering of the six independent samples was obtained for selected groups of genes using GeneSpring 7.3. Each colored box 

represents the normalized expression level of a given gene in each sample and is colored according to the fold change. B. Line Graph. 

Represents a two-fold differential gene expression between the AMC and RMC. The lines in red represent gens upregulated and those in blue 

represent genes downregulated in AMC in comparison to RMC. C. Volcano Plot. Within the lateral quadrants (red and blue box) are the genes 

with two-fold difference and P Value < 0.05. These genes were chosen for generation of functional group lists. D.Validation of Microarray 

Profile. Histogram shows the qRTPCR validation of two AMC-specific genes (Runxltl and Sept9) and two RMC-specific genes (Sept4 and Mbp). 
k J 



everse: 5 -GAGCTGATGCAGGGAAG-3 Mbp are for- 
ward: 5 -TACTTGGCCACAGCAAGTACC-3 ' and re- 
verse: 5 -GGGTGTACGAGGTGTCACAAT-3 Gapdh are 
forward: 5 -ACATGCCGCCTGGAGAAACCTGCCA-3 ' 



and reverse: 5-TGCCAGCCCCAGCATCAAAGGTGGA- 
3'. The primer sequences used for the data reported in the 
supplementary figure are listed in Additional file 1: Sheet 
SI. After pre-incubation at 95°C for 15 min, the polymerase 
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Table 1 Top 25 highly-expressed genes in AMC based on fold change. Both MATLAB and GeneSpring analysis are 
represented 



Genes from MATLAB analysis 



Gene Symbol Gene Title 



Sla 

Dpysl3 
Sox4 



Rprm 

Satb2 
Dcx 

Crmpl 
Syt16 
Cct8 



Src-like adaptor 



Fold Function 
Change 

73.34 protein binding 



dihydropyrimidinase-like 3 31.01 nervous system 

development 



SRY (sex determining 
region Y)-box 4 



reprimo, TP53 dependent 
G2 arrest mediator 
candidate 

SATB homeobox 2 



Doublecortin 



collapsin response 
mediator protein 1 

synaptotagmin XVI 



chaperonin containing 
Tcpl, subunit 8 (theta) 

RGD1310352 similar to HTGN29 

protein; keratinocytes 
associated 

transmembrane protein 2 



Appbp2 
Cxcr4 

Ect2 

Hmgb3 
Hs3st5 

Cfl2 
Tnc 

Napl 13 
Mex3b 
Smarcal 



28.33 pro-B cell 

differentiation 



19.13 cell cycle arrest 

18.05 negative regulation 
of transcription 

16.64 neuronmigration 

15.50 neuron development 

13.53 protein binding 

13.53 protein folding 

13.42 Unknown 

13.28 Transport 



receptor 4 

epithelial cell transforming 12.50 
sequence 2 oncogene 

high mobility group box 3 12.47 

heparan sulfate 12.32 
(glucosamine) 
3-O-sulfotransferase 5 



cofilin 2, muscle 
Tenascin C 



nucleosome assembly 
protein 1 -like 3 

mex3 homolog B 
(C. elegans) 

SWI/SNF related, matrix 
associated, actin dependent 
regulator of chromatin, 
subfamily a, member 1 



12.21 
11.93 

10.88 
10.78 



Genes from Genespring 

Gene Symbol Gene Title 

LRRGT00193 Unknown 

Rnfl 52 ring finger protein 152 



Syncrip Synaptotagmin 
binding, 

cytoplasmic RNA 
interacting protein 

Chst8 carbohydrate (N- 

acetylgalactosamine 4-0) 
sulfotransferase 8 

Rpl28 ribosomal protein L28 



Slc16a7 solute carrier family 16, 

member 7 

(monocarboxylic acid 
transporter 2) 

Crmpl collapsin response 

mediator protein 1 

LOC1 001 58225 hypothetical protein 
LOC1 001 58225 



analysis 

Fold 
Change 

574.79 

35.04 

26.62 



amyloid beta precursor 
protein (cytoplasmic tail) 
binding protein 2 

chemokine (C-X-C motif) 12.76 ameboidal cell 



Bmp2 



Dcx 



Mab21M 



Hmgb3 



migration 



cell morphogenesis Bazla 



negative regulation 
ofB cell differentiation 

protein amino 
acid sulfation 



protein binding 



negative regulation 
of cell adhesion 



nucleosome 
assembly 

RNA binding 



10.72 brain development 



Sox4 

LOC688455 

Cxcr4 
Spint2 

Zfml 

Nradd 

Bcl7c 



bone morphogenetic 
protein 2 

doublecortin 



mab-21-like 1 
(C. elegans) 

high mobility group 
box 3 



bromodomain 
adjacent to zinc finger 
domain, 1A 

SRY (sex determining 
region Y)-box 4 

hypothetical protein 
LOC688455 



26.17 

19.32 
18.85 

18.54 
17.44 
16.43 
15.26 



13.43 



chemokine (C-X-C motif) 1 2.92 
receptor 4 

serine peptidase 12.49 
inhibitor, Kunitz type, 2 

zinc finger, matrin-like 12.20 



neurotrophin receptor 11.65 
associated death domain 



Function 

Unknown 
Unknown 

mRNA 
processing 

sulfur 

metabolic 

process 

translation 
transport 



neuron 
development 

Unknown 

ossification 

neuron 
migration 



14.54 positive 

regulation of 
cell proliferation 

14.14 negative 

regulation of 
myeloid cell 
differentiation 

13.89 protein binding 



pro-B cell 
differentiation 



12.93 Unknown 



ameboidal 
cell migration 

serine-type 
endopeptidase 
inhibitor activity 

nucleic acid 
binding 

signal 

transduction 



B-cell CLL/lymphoma 7 C 1 1.44 Unknown 
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Table 1 Top 25 highly-expressed genes in AMC based on fold change. Both MATLAB and GeneSpring analysis are 
represented (Continued) 



Sh3bgrl 


SH3 domain binding 
glutamic acid-rich 
protein like 


9.74 


Unknown 


Tmeffl 


transmembrane protein 
with EGF-like and two 
fol I istati n-li ke domains 1 


11.43 


multicellular 

organismal 

development 


Ascc3l1 


Activating signal 
cointegrator 1 
complex subunit 3-like 1 


9.58 


RNA splicing 


Mtusl 


mitochondrial tumor 
suppressor 1 


11.12 


cell cycle 


Bex4 


brain expressed gene 4 


9.50 


Nucleus 


Adfp 


Adipose differentiation 
related protein 


10.85 


response to 
organic cyclic 
substance 


Fam164a 


family with sequence 
similarity 164, member A 


9.49 


Unknown 


Ttk 


Ttk protein kinase 


10.65 


protein amino 
acid fig7 
phosphorylation 


LOC294446 


similar to Myristoylated 
alanine-rich C-kinase 
substrate (MARCKS) 
(ACAMP-81) 


9.44 


calmodulin binding 


Lrrc20 


leucine rich repeat 
containing 20 


10.63 


protein binding 



chain reaction (PCR) was performed as follows: 45 cycles of 
denaturation at 94°C for 15 s, annealing at 57°C for 25 s, 
and elongation at 72°C for 15 s. 

Results 

Laser-capture microdissection of microglial cells from 
the corpus callosum of 5-day and 4-week old rat brain. 

To compare the gene expression profiles of AMC and 
RMC, we stained both microglial cell types with 
peroxidase-conjugated lectin and isolated them from the 
CC of 5-day and 4-week old rat brain respectively. LCM of 
AMC and RMC from the CC of 5 -day old rat brain has 
been shown in Figure 1A-F. Lectin staining has been 
widely used to selectively stain microglia for study of 
microglial development in the CNS [29-31]. The cells iso- 
lated by LCM were further confirmed to be microglia 
since the mRNA expression of oligodendrocyte (CNPase), 
astrocyte (GFAP) and endothelial cell-specific genes 
(Vimentin) was undetectable (data not shown). 

cDNA microarray and generation of gene lists specific to 
AMC and RMC 

To identify the genes that are differentially expressed be- 
tween AMC and RMC, we extracted total RNA from AMC 
and RMC and carried out cDNA microarray using Rat 
Genome 230 2.0 array (Affymetrix). Each sample contained 
RNA from six hundred laser-captured microglial cells. To 
ensure gene expression consistency between samples within 
the groups, we determined the Pearsons rank correlation 
coefficient after normalizing the raw expression data 
(Figure 1G). The gene expression profile from the samples 
of same group showed a very high correlation of 0.97 ± 0.03 
while, a relatively lower correlation value of 0.87 ± 0.03 was 
observed between samples of different groups. A high cor- 
relation coefficient of above 0.8 between the AMC and 
RMC may be due to the fact that the comparison is 



between the gene expression profiles of the same cell type., 
i.e. microglia regardless of the differences in age (5-day and 
4-week rat brain) and morphology (amoeboid and 
ramified). 

Agglomerative average-linked hierarchical clustering 
was performed and genes showing over two-fold differ- 
ential expression between AMC and RMC were identi- 
fied using GeneSpring 7.3 (Figure 2A & B). About 800 
genes were found to be differentially regulated in the 
two sample groups - 537 with upregulation and 258 with 
downregulation in AMC. A high number of differentially 
expressed genes identified were either novel or did not 
have any functional annotation. In view of this, a list of 
genes with known functional annotations was generated 
using the statistical functions in the Bioinformatics Tool- 
box in MATLAB R2009a [32]. This list was generated 
using a less stringent filtering (p value < 0.05, in contrast 
to the GeneSpring list which has a p value < 0.01) and 
contained close to 1400 upregulated genes and 700 
downregulated genes in AMC compared to RMC 
(Figure 2C). 

Functional categorization of AMC and RMC 

The genes with highest fold change values in AMC and 
RMC (Table 1 and 2) clearly delineate their functions and 
residing environment. For example, the AMC express 
genes that are shown to be involved in nervous system de- 
velopment (such as Dpysl3, Crmpl and Smarcal) [33,34], 
immune system development (such as Hmgb3 and Sla) 
[35-38], cell migration during neurodevelopment (Dcx) 
[39] and the immune response as well in migration of 
microglia (Cxcr4, a chemokine receptor) [40-42]. The 
finding of expression of some genes that are known to be 
neuron-specific (such as Dcx), is interesting and has been 
further confirmed by immunohistochemical analysis 
which revealed the expression of Dcx by the AMC in the 
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Table 2 Top 25 genes in RMC based on fold change. Both MATLAB and GeneSpring analysis are represented 



Genes from MATLAB analysis Genes from Genespring analysis 



Gene 
Symbol 


Gene Title 


Fold 
Change 


Function 


Gene 
Symbol 


Gene Title 


Fold 
Change 


Function 


Mobp 


myelin-associated 
oligodendrocyte 
basic protein 


623.18 


nervous system 
development 


Mog 


myelin oligodendrocyte 
glycoprotein 


26.34 


cell adhesion 


Mog 


myelin oligodendrocyte 
glycoprotein 


39.14 


cell adhesion 


Slc5a1 1 


solute carrier family 5 
(sodium/glucose 
cotransporter), 
member 1 1 


24.41 


antigen processing 
and presentation 


Mbp 


myelin basic protein 


21.57 


myelination 


Mbp 


myelin basic protein 


22.86 


myelination 


Robo3 


roundabout homolog 3 
(Drosophila) 


1 7.39 


neuron migration 


Cyp3a9 


cytochrome P450, 
family 3, subfamily a, 
polypeptide 9 


1 9.31 


sensory perception 
of smell 


Mai 


mal, T-celldifferentiation 
protein 


16.13 


intracellular protein 
transport 


Opalin 


oligodendrocytic myelin 
paranodal and 
inner loop protein 


15.59 


Golgi apparatus 


Dnah12 


dynein, axonemal, 
heavy polypeptide 12 


12.77 


microtubule-based 
movement 


Plpl 


proteolipid protein 1 


12.00 


glial cell 
differentiation 


Slc5a11 


solute carrier family 5 
(sodium/glucose 
cotransporter), member 1 ' 


12.34 

I 


antigen processing 
and presentation 


Camk2a 


calcium/calmodulin- 
dependent protein kinase 
II alpha 


9.59 


G1/S transition 
of mitotic cell 
cycle 


Hapln2 


hyaluronan and 
proteoglycan link 
protein 2 


9.84 


cell adhesion 


Lgi4 


leucine-rich repeat 
LGI family, member 4 


9.39 


neuron maturation 


Plp1 


proteolipid 
protein 1 


9.00 


glial cell 
differentiation 


Cpne9 


copine family 
member IX 


8.71 


unknown 


Ermn 


ermin, ERM-like protein 


8.82 


morphogenesis of a 
branching structure 


Znf76 


zinc finger protein 76 
(expressed in testis) 


8.55 


transcription 


Tnnc2 


troponin C type 2 (fast) 


8.60 


skeletal muscle 
contraction 


Rhpnl 


rhophilin, Rho GTPase 
binding protein 1 


8.35 


signal transduction 


Mag 


myelin-associated 
glycoprotein 


8.12 


cell adhesion 


Gng13 


guanine nucleotide 8.16 
binding protein (G protein), 
gamma 13 


G-protein coupled 
receptor protein 
signaling pathway 


Lgi4 


leucine-rich repeat 
LGI family, member 4 


8.10 


neuron maturation 


Mag 


myelin-associated 
glycoprotein 


7.99 


cell adhesion 


Extll 


exostoses (multiple)-like 1 


8.09 


protein binding 


Herc6 


hect domain and RLD 6 


7.71 


protein modification 
process 


Aldh3b1 


aldehyde 
dehydrogenase 3 
family, member B1 


7.81 


cellular aldehyde 
metabolic process 


Chn1 


Chimerin (chimaerin) 1 


6.31 


signal transduction 


Hhatl 


hedgehog 
acyltransferase-like 


7.79 


negative regulation 
of N-terminal protein 
palmitoylation 


Sept4 


septin 4 


5.81 


cell cycle 


Gng13 


guanine nucleotide 

binding protein 

(G protein), gamma 13 


7.44 


G-protein coupled 
receptor protein 
signaling pathway 


Ccdc37 


coiled-coil domain 
containing 37 


5.66 


unknown 


Akap8l 


A kinase (PRKA) anchor 
protein 8-like 


7.30 


DNA binding 


Cmtm5 


CKLF-like MARVEL 
transmembrane domain 
containing 5 


5.30 


membrane 


Cpne9 


copine family 
member IX 


7.24 


Unknown 


Amhr2 


anti-Mullerian 
hormone receptor, 
type II 


5.08 


protein amino 

acid phosphorylation 


Carl 2 


Carbonic anyhydrase 12 


7.21 


one-carbon 
metabolic process 


Sptbn4 


spectrin, beta, 
non-erythrocytic 4 


4.84 


axonogenesis 


Srpk3 


SFRS protein kinase 3 


7.11 


protein amino acid 
phosphorylation 


Gsta3 


glutathione 
S-transferase A3 


4.55 


glutathione 
metabolic process 
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Table 2 Top 25 genes in RMC based on fold change. Both MATLAB and GeneSpring analysis are represented 

(Continued) 



Casd 

S1pr5 
Gpd1 

Itsn2 



cancer susceptibility 6.90 
candidate 1 



sphingosine-1 -phosphate 6.82 
receptor 5 

glycerol-3-phosphate 6.75 

dehydrogenase 1 

(soluble) 



intersectin 2 



6.67 



Unknown 



signal transduction 

carbohydrate 
metabolic process 



regulation of Rho 
protein signal 
transduction 



Grin2c glutamate receptor, 4.53 
ionotropic, N-methyl 
D-aspartate 2 C 

Fah fumarylacetoacetate 4.50 

hydrolase 

Ggn gametogenetin 4.49 



Ntsr2 neurotensin 4.13 
receptor 2 



startle response 



arginine catabolic 
process 

multicellular 

organismal 

development 

signal transduction 
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Figure 3 A-F. Differential immunoexpression of Dcx andRunxltl in AMC and RMC. Confocal images showing the immunoexpression of 
Dcx (B; red) and its co-localization (C) in 0X42 (A; green) labeled AMC. Immunoexpression of Runxltl (E, H; red) and its co-localization (F, I) in 
0X42 (D, G; green) labeled AMC and RMC in the CC from 5-day (5D) and 4-week (4 W) old rat brain was also observed. Runxltl 
immunoexpression is undetectable in RMC (I) compared to that in the AMC (F). (DAPI - blue). Scale bars: A-C 50 urn, D-l 10 urn. 
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Figure 4 Functional clusters to highly-expressed AMC and RMC genes. Heat map shows the top 25 AMC and RMC (arranged according to 
fold change) and their involvement in major cellular functions. Red shading indicates AMC genes and green shading indicates RMC genes. 



CC of 5 day old rat brain (Figure 3A-C). In addition, 
mRNA expression of Dcx was detected in the BV-2 micro- 
glia (a murine microglia cell line) by RTPCR (Additional file 
2: Figure SI). 

Interestingly, the RMC express genes involved in myelin- 
ation (such as Mbp). Mbp-like proteins, also known as 
Golli proteins have previously been shown to be localized 
in human microglia at 22 weeks postnatally [43]. MBP 
mRNA was found to be expressed by laser captured-RMC 
(Figure 2D) and BV-2 microglia (Additional file 2: Figure 
SI). Other myelin-related genes like Plpl [44] and Lgi4 [45] 
have also been found with high expression values in the 
RMC. Plpl and Lgi4 were found to be expressed by non- 
myelinating cells such as the Bergmann glia in cerebellum 
of the developing mouse brain [46]. 

On sorting the genes based on p values (Additional file 3: 
Sheet S2), we found several genes that are specific to AMC 
such as, genes involved in transcriptional repression (Mbdl 
which binds to methylated sites on DNA) [47,48], vesicular 
trafficking (Snx6, a component of the retromer complex) 
[49], and microtubule depolymerization (Stmnl) [50]. RMC 
express genes involved in immune functions such as RT1- 



A2, which is the MHC of rat and Clql3, a protein of the 
complement system [51,52], calcium ion signaling pathway 
protein, Camk2 [53,54] and sodium dependent glucose 
transporter gene Slc5all, known to interact with immune- 
related genes [55]. 

bone-marrow derived circulating monocytes [2]. The 
availability of gene expression profiles for circulating 
monocytes prompted us to compare and study the func- 
tional similarities of AMC and RMC to that of circulat- 
ing monocytes [28] (Refer Table 3 & 4 for pathways, 
Additional file 7: Sheet S6 for full gene list). Monocytic 
genes highly expressed by AMC (Cxcr4, Csk and Racl) 
are known to be involved in several disease pathways 
such as Parkinsons disease, Huntington's disease and 
HIV infection, phagocytosis and chemokine signaling 
pathways. Further, the monocytic genes (Hla-c, Cd74, 
Cd302, Lspl and Runx3) expressed by RMC are involved 
in antigen presentation and lysosome-related functions. 
The peripheral monocytic genes expressed by the AMC 
and RMC may be further investigated if they reflect the 
microglia property in pathological brain. These genes 
may also provide deeper insights into the similarities 
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Figure 5 A. Functional pathways. Graph shows the number of AMC and RMC genes enriched in different functions. It can be noted that the 
AMC apart from proliferation and differentiation express a high number of genes involved in cell death whereas the RMC express a high number 
of cytoskeletal genes. B. Sternness of AMC and RMC. Venn diagrams represent the share of Embryonic (ESC), Neural (NSC) and Hematopoietic 
(HSC) stem cell genes in AMC and RMC. The AMC express a high number of stem cell genes compared to the RMC. The RMC express more HSC- 
specific genes than NSC- and ESC- specific genes. 



between peripheral monocytes and microglia and the re- 
appearance of these genes during infection or neurode- 
generation in the activated microglia may be critical for 
immune response. 

Expression of neural, embryonic and hematopoietic stem 
cell specific genes by AMC and RMC 

Both AMC and RMC express a number of stem cell 
specific genes [27] (Table 5 &6, Additional file 8: 
Sheet S7). However, the RMC express a lesser num- 
ber of neural stem cell (NSC) and embryonic stem 
cell (ESC) specific genes, compared with the AMC. In 
spite of this, higher percentage of hematopoietic stem 
cell specific genes (HSC) was found to be expressed 
in RMC (Figure 5B). Our finding that microglia retain 



HSC specific properties even in the adult brain is 
suggestive of their hematopoietic lineage. 

Pathway analysis and validation of differentially 
expressed genes in AMC and RMC 

The pathway analysis (Additional file 9: Figure S2) revealed 
novel molecular networks involving several signaling mole- 
cules and pathways within microglia. In order to validate 
the results obtained from pathway analysis, we randomly 
selected three transcription factors which are highly 
expressed in the AMC. They are: Sox4 and Soxll which 
are SRY-related HMG-box family of transcription factors 
[57] and Runt-related transcription factor 1; translocated 
to, 1 (cyclinD-related) (Runxltl), a member of the ETO 
gene family of transcriptional co-repressors [58]. 
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Table 3 Pathways involving monocytic genes expressed by AMC 



Pathway 


Nr. of 


P-Value 


Benjamini 




Genes 




corrected P-value 


REACT_17015:Metabolism of proteins 


37 


2.00E-16 


1.20E-14 


REACT_71:Gene Expression 


42 


3.30E-13 


8.60E-12 


REACT_1 762:3' -UTR- mediated translational regulation 


24 


6.20E-13 


1.10E-11 


Ribosome 


17 


1.30E-10 


1.40E-08 


REACT_6167:lnfluenza Infection 


24 


1 .90E-09 


2.50E-08 


Proteasome 


10 


1 .20E-06 


6.90E-05 


Pathogenic Escherichia coli infection 


10 


6.70E-06 


2.50E-04 


P00029:Huntington disease 


16 


7.20E-06 


5.80E-04 


REACT_13635:Regulation of activated PAK-2p34 by proteasome mediated degradation 


10 


1 .50E-05 


1.50E-04 


P00049:Parkinson disease 


12 


3.40E-05 


1.40E-03 


REACT_9035:APC/C:Cdh1 -mediated degradation of Skp2 


11 


3.80E-05 


3.30E-04 


REACT_6850:Cdc20:Phospho-APC/C mediated degradation of Cyclin A 


11 


7.60E-05 


5.70E-04 


REACT_6185:HIV Infection 


18 


1 .70E-04 


1.10E-03 


REACTJ 1045:Signaling by Wnt 


10 


2.20E-04 


1.30E-03 


P00016:Cytoskeletal regulation by Rho GTPase 


10 


1 .40E-03 


3.70E-02 


Fc gamma R-mediated phagocytosis 


9 


1 .90E-03 


5.10E-02 


P00018:EGF receptor signaling pathway 


11 


2.40E-03 


4.80E-02 


P00034:lntegrin signalling pathway 


14 


2.70E-03 


4.20E-02 


Neurotrophin signaling pathway 


10 


2.70E-03 


6.00E-02 


P00060:Ubiquitin proteasome pathway 


8 


4.30E-03 


5.70E-02 


REACT_578:Apoptosis 


12 


4.90E-03 


2.50E-02 


Chemokine signaling pathway 


12 


4.90E-03 


8.80E-02 


REACTJ 538:Cell Cycle Checkpoints 


11 


5.30E-03 


2.50E-02 


P00021:FGF signaling pathway 


10 


5.80E-03 


6.50E-02 


REACT_383:DNA Replication 


10 


6.20E-03 


2.70E-02 



Runxltl, by forming a fusion protein with Runxl, an- 
other member of RUNX family, leads to self-renewal of 
human monocytic cells thereby impairing differentiation of 
these cells [59]. Certain genes known to be downregulated 
by the Runxl -Runxltl transcription factor complex such 
as Socsl, Csfl, and Runx3 are highly expressed by RMC 
[60] (Additional file 5: Sheet S4). Further, a transcriptional 
dysregulation caused by this fusion protein was found to 
cause the over-expression of Sox4 in human progenitor 
cells [60]. In an earlier study, Sox4 deficient mice exhibited 
proliferation-defective pro-B cells [61]. Similar networks 
might function in the AMC and therefore warrants further 
investigation. Our immunohistochemical analysis revealed 
that Sox4 (Figure 6A-C) is highly expressed in the nucleus 
and cytoplasm of AMC. Similarly, Soxll (Figure 6D-I) and 
Runxltl (Figure 3D-I) are expressed in the AMC, but 
hardly detectable in the RMC thus validating our micro- 
array results. Further, quantitative real time RTPCR for 
Runxltl using RNA extracted from LCM-captured AMC 
and RMC showed a very high expression of Runxltl in 
AMC when compared to RMC (Figure 2D). 



Differential expression of septin genes in AMC and RMC 

Septins are a family of multifunctional proteins involved in 
cytoskeletal organisation and cell division [62-64]. They 
have also been implicated in tumorigenesis and neurode- 
generation [65]. In the present study, Septin genes were 
found to be differentially expressed in AMC and RMC. 



Table 4 Pathways involving monocytic genes expressed 
by RMC 



Pathway 


Nr. of 
Genes 


P-Value 


Benjamini 
corrected 
P-value 


Lysosome 


5 


2.00E-03 


1.00E-01 


Antigen processing 
and presentation 


4 


6.60E-03 


1.60E-01 


70.Signal_peptides_(MHC)_class_l 
_molecules 


2 


3.30E-02 


1.60E-01 


Other glycan degradation 


2 


7.30E-02 


7.40E-01 


1 14.Genomic_reformatting_Brain 
Jschemia 


2 


9.70E-02 


2.30E-01 
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Table 5 Stem-cell genes enriched in AMC 


Gene Symbol 


Gene Title 


Fold Changes 


Function 


Rprm 


reprimo, TP53 dependent G2 arrest mediator candidate 


1 9.1 25 


cell cycle arrest 


Cct8 


chaperonin containing Tcp1, subunit 8 (theta) 


1 3.525 


protein folding 


nbU I o I Uojz 


similar to HTGN29 protein; keratinocytes associated 
transmembrane protein 2 


1 3.4 1 y 


Unknown 


Ect2 


epithelial cell transforming sequence 2 oncogene 


12.498 


cell morphogenesis 


Cf\2 


cofilin 2, muscle 


12.211 


protein binding 




SH3 domain binding glutamic acid-rich protein like 


9.7435 


Unknown 


Bex4 


brain expressed gene 4 


9.4991 


Unknown 


Cfl1 


cofilin 1, non-muscle 


9.3235 


cytokinesis 


Cd24 


CD24 molecule 


9.2544 


response to hypoxia 


Ns5atp9 


NS5A (hepatitis C virus) transactivated protein 9 


9.1968 


Unknown 


Hnrnpab 


heterogeneous nuclear ribonucleoprotein A/B 


8.8712 


epithelial to mesenchymal transition 


Mycn 


v-mycmyelocytomatosis viral related oncogene, 
neuroblastoma derived (avian) 


8.6884 


regulation of transcription, 
DNA-dependent 


LOC294446 


similar to Myristoylated alanine-rich C-kinase substrate 
(MARCKS) (ACAMP-81) 


8.5949 


actin binding 


Pex13 


peroxisomal biogenesis factor 13 


8.4726 


fatty acid alpha-oxidation 


Maoa 


monoamine oxidase A 


7.8728 


catecholamine metabolic process 


Tmeffl 


transmembrane protein with EGF-like and 
two follistatin-like domains 1 


7.7994 


multicellular organismal development 


Ube2e3 


ubiquitin-conjugating enzyme E2E 3, UBC4/5 homolog (yeast) 


7.5905 


modification-dependent protein 
catabolic process 


Rab10 


RAB10, member RAS oncogene family 


7.5637 


regulation of transcription, 

PlKIA rlononrlont 




r\Ur 1 lUUbyidllUI 1 IdLlUl 


/ .DDZZ 


ll dl IbpUl I 


Maprel 


microtubule-associated protein, RP/EB family, member 1 


7.0496 


cell cycle 


Psme3 


proteasome (prosome, macropain) activator subunit 3 


7.0045 


cell adhesion 


Tmem43 


transmembrane protein 43 


6.9166 


Unknown 


March5 


membrane-associated ring finger (C3HC4) 5 


6.9148 


zinc ion binding 


Sqle 


Squaleneepoxidase 


6.7699 


cellular aromatic compound metabolic 
process 


Apafl 


apoptotic peptidase activating factor 1 


6.6501 


neural tube closure 



Sept3, 6, 9 and 11 were expressed in AMC whereas Sept4 
and 8 expressed in RMG The expression of some of these 
genes was further confirmed by immunohistochemical and 
quantitative real time RTPCR analysis which showed Sept9 
immunoexpression and mRNA expression in the AMC 
(Figure 7A-F and Figure 2D) and Sept4 expression in the 
RMC (Figure 7G-L and Figure 2D). Since the role of Sep- 
tins in microglial functioning has not yet been investigated, 
studies on the Septin family may further our knowledge on 
the cytoskeletal dynamics involved in proliferation, migra- 
tion and activation of AMC and RMC. 

Microglia in response to trauma or neurodegenerative 
stimuli exhibit upregulation of proinflammatory cyto- 
kines and chemokines [66]. In the present study, both 
AMC and RMC exhibited relatively low expression in- 
tensities for most of the cytokines such as TNF-a and 



interleukins and chemokines such as Cxcl3, Cxcll2 and 
Ccl2 (Additional file 10: Sheet S8). 

Discussion 

Microglia in the healthy brain express low levels of 
cytokines and chemokines 

This study is a novel attempt to examine the global gene 
expression profile of microglia in situ and to functionally 
distinguish the two distinct microglial phenotypes, namely, 
AMC and RMC. A noteworthy feature of this transcrip- 
tome profile was that the expression of cytokines and che- 
mokines in both AMC and RMC was hardly detectable 
which is in agreement with previous studies [24]. It has 
been widely shown that the untreated microglia in culture 
produce some amount of proinflammatory cytokines and 
chemokines, indicating that culture media stimulate the 
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Table 6 Stem-cell genes enriched in RMC 



Gene Symbol 


Gene Title 


Fold Changes 


Function 


RGD1 307882 


similar to CG9346-PA 


2.9758 


RNA processing 


Mill 


myeloid/lymphoid or mixed-lineage leukemia 1 


2.9906 


DNA repair 


Fermt3 


fermitin family homolog 3 (Drosophila) 


3.0629 


protein binding 


Ppfibp2 


PTPRF interacting protein, binding protein 2 (liprin beta 2) 


3.2355 


DNA integration 


Mtfl 


metal-regulatory transcription factor 1 


3.2416 


regulation of transcription, 
DNA-dependent 


Egr1 


early growth response 1 


3.2473 


negative regulation of transcription 
from RNA polymerase II promoter 


Map7 


microtubule-associated protein 7 


3.263 


cell morphogenesis 


Sfrsl 1 


splicing factor, arginine/serine-rich 1 1 


3.3325 


nuclear mRNA splicing, via 
spliceosome 


LOoUz 


f~VYMY} r>nr\\r\/~i \\r\ 

LUoUz molecule 


Q 37/1 ~) 


binding 


C tnnal 1 
L.LI II Idl I 


LdLclllll ^LdUllcllll dSSUUdlcU piOlclNJ, dipild II Ktr I 


3 A~\ 11 
Orr \ A/ 


unknown 




r\/tr>rhrr>mQ D/l QP\ A PA 

cy loci II oi t ie r^ou *\ ro 




leukotriene metabolic process 


Mlec 


Malectin 


3.5937 


carbohydrate metabolic process 


Smarca2 


SWI/SNF related, matrix associated, actin 

Ucptrl lUfcrl I L IcLjUldLUI Ul LI IIUI 1 IdLII I, bUUIdl 1 Illy d, IllclllUcI Z 


3.663 


negative regulation of cell 

piUlllcldUUl I 


OCI I I U 


sodium channel, voltage-gated, type 1, beta 




transport 


I 1 1 lUll 1 1 1 


Lldl Iblllcl IlUldl lc DMA INNIUILUI IIIULII LUI lldll III iy I 


Q 7QQO 


unknown 


r CKz 


phosphoenolpyruvatecarboxykinase 2 (mitochondrial) 




gluconeogenesis 


Tnni3 


LfUpUNIN I Lypt: D ^LdlUldLJ 


A f)f)AR 


VdbLUIULjcl Icblb 


Mpfh 

1 NCI 1 1 


I ItfUlUllldl I ltrl IL, NtrdVy pUiypcpilUc 


a ?n^^ 


1 1 IILI ULUUUIc Ly ICJSKclclON 

organization 


Prpf38b 


PRP38 pre-mRNA processing factor 38 (yeast) domain containing B 


4.3721 


mRNA processing 


Tnntl 


troponin T type 1 (skeletal, slow) 


4.4062 


skeletal muscle contraction 


Zranb2 


zinc finger, RAN-binding domain containing 2 


4.598 


mRNA processing 


Zrsrl 


zinc finger (CCCH type), RNA binding motif and serine/arginine rich 1 


5.0348 


nucleotide binding 


Lrrc23 


leucine rich repeat containing 23 


5.1395 


protein binding 


Akap8l 


A kinase (PRKA) anchor protein 8-like 


5.6932 


DNA binding 


Tnnc2 


troponin C type 2 (fast) 


8.6024 


skeletal muscle contraction 



microglial cells. Significantly, the low expression of cyto- 
kines and chemokines in both types of microglia in the 
present study appears to mimic the transcriptome status of 
normal microglia in healthy brain in vivo. 

AMC express genes involved in cell cycle process and 
migration whereas RMC express genes involved in 
synaptic integrity and neuronal maturation 

AMC from the first week of postnatal rat brains have a high 
proliferative capacity [2]. During development, about two- 
thirds of AMC undergo apoptosis and the rest transform 
into RMC [67,68]. In accordance with this, microarray ana- 
lysis in the present study revealed a high expression of cell 
proliferation/cell cycle-related genes such as Myc and 
CyclinA2, CyclinB2 and CyclinDl (Additional file 6: Sheet 
S5) [69,70] and genes involved in cell death (Figure 5A & 
Additional file 6: Sheet S5) namely, Casp2,Casp3 [71,72] 
and Apafl (Glial apoptosis related gene) in the AMC [73]. It 



is striking that the AMC express Dcx, a protein known to 
be a marker for migrating neurons [74]. It may be worth in- 
vestigating the role of Dcx in migration of AMC in the 
early postnatal brain. 

On the other hand, the RMC, apart from cell homeosta- 
sis and glial development, appear to contribute to synaptic 
transmission as they express genes such as Grin2c, SlOOb 
and Camk2a (Additional file 5: Sheet S4) [75-77]. This is 
interesting and supports the recent experimental studies 
showing the role of microglia in the maintenance and 
modifications of synaptic integrity in the healthy brain 
[78,79]. Further, Grin2c, a subunit of NMDA receptor 
complex is expressed by the microglial cells in the CC and 
has been shown to be functionally important in microglia- 
mediated neuroinflammation [80,81]. SlOOb, a calcium ion 
binding protein, is also expressed by microglia and relo- 
cates around phagosomes during microglial activation and 
phagocytosis [82,83]. 
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Figure 6 A-l. Differential immunoexpression of Sox4 and Soxll in AMC and RMC. Confocal images showing a high immunoexpression of 
Sox4 (B; red) and its co-localization (C) in 0X42 (A; green) labeled AMC. Immunoexpression of Soxll (E, H; red) and its co-localization (F, I) in 
0X42 (D, G; green) labeled AMC and RMC in the CC from 5-day (5D) and 4-week (4 W) old rat brain was also observed. Sox1 limmuno expression 
is undetectable in RMC (I) compared to that in the AMC (F). (DAPI - blue).Scale bars: A-l 10 urn. 



RMC express myelin basic protein (Mbp) which 
encodes two families of proteins i.e., classic Mbps and 
Golli-Mbps. Classic Mbps serve as the major protein 
constituent of myelin in the central and peripheral ner- 
vous system whereas, the Golli proteins are known to be 
broadly distributed in the brain, particularly in the 
microglia during normal development and inflammation 
[84] and involved in the interaction between microglia 
and oligodendrocyte precursor cells during multiple 
sclerosis [85]. Taken together, these results reveal that 
microglia are not only involved in immune response and 
phagocytosis but also play diverse roles in healthy brain. 

Both AMC and RMC express cytoskeleton-related genes 

Regulation of cytoskeletal dynamics is important to both 
microglial migration and ramification [86-88]. Apart from 



cytoskeletal structural proteins such as tubulins and actin, 
we found that the AMC express cytoskeleton-associated 
Crmp family proteins [89] such as Crmp-1, Dpysl3 and 
Dpysl5 and Septin family proteins such as Sept9 and 11 
[90,91]. Sep tins are implicated in cytoskeletal processes 
such as vesicular trafficking [92]. These cytoskeleton- 
associated proteins may therefore explain the migration 
and phagocytosis of AMC during normal development 
and pathology. In the present study, AMC express Sept9 
but not Sept4 whereas, RMC express Sept4, but not Sept9, 
indicating differential roles of Septin family genes in AMC 
and RMC. Sept4 has been recently shown to be involved 
in cortical neuron migration [93]. Absence of Sept4 
immunoexpression in the AMC and its high expression in 
the RMC is suggestive of an important role for this protein 
in microglial transformation during development. 
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Figure 7 A-L. Immunoexpression of Sept9 and Sept4 in RMC. Confocal images showing the immunoexpression of Sept9 and Sept4 (B, E, H, 
K; red) and their co-localization (C, F, I, L) in 0X42 (A, D, G, J; green) labeled AMC and RMC in the CC from 5-day (5D) and 4-week (4 W) old rat 
brain.Sept9 immunoexpression is undetectable in RMC (F) compared to that in the AMC (C) whereas Sept4 is undetectable in AMC (I) compared 
to that in the RMC (L). (DAPI - blue).Scale bars: A-L 50 urn. 



Expression of monocyte- and stem cell-specific genes by 
AMC and RMC indicates their sternness and origin 

Recent studies have proven that microglia originate from 
the mesenchymal progenitor cells at the yolk sac [56]. 
However, microarray studies of various hematopoetic and 
non-hematopoetic cell types revealed a close relationship 
between the gene expression profiles of microglia and 
bone- marrow derived macrophages [94] which are known 



to differentiate from circulating monocytes. Therefore, we 
sought to identify the monocyte-specific genes expressed 
by AMC and RMC. AMC express several monocyte- 
specific genes including Mcll and Id2. Mcll is associated 
with cell viability and differentiation of myeloid cells 
which include monocytes and macrophages [95,96] and 
Id2, a negative regulator of basic helix loop helix tran- 
scription factors, is involved in the differentiation of 
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myeloid cells [97]. A recent study demonstrated that Id2 is 
required for bone morphogenic protein (Bmp) -mediated 
differentiation of microglia into Map2 + neurons and Gfap + 
astrocytes [98] suggesting that this gene may promote 
microglial trans -differentiation. Both Mcll and Id2 have 
been shown to be involved in cell differentiation and their 
high expression in AMC explains the role of these genes in 
promoting the maturation of AMC and its transformation 
into RMC 

On the other hand, RMC exhibited increased expres- 
sion of Lspl, which binds to the cytoskeleton and is 
known to be a marker for leucocytes [99]. Overexpres- 
sion of Lspl in neutrophils was associated to defective 
actin polymerization which render these cells immotile 
[100]. Further, overexpression of Lspl in a highly motile 
melanoma cell line led to formation of hair-like projec- 
tions. Thus, upregulation of Lspl in the RMC, compared 
to the AMC may explain the role of this gene in motility 
and ramification of RMC [101] which are the resident 
population in the adult brain parenchyma. 

Both AMC and RMC express stem cell-specific (ESC, 
NSC and HSC) genes, indicating their sternness and sug- 
gesting that microglia may undergo trans -differentiation. 
The RMC expressed a high percentage of HSC specific 
genes in comparison to ESC and NSC specific genes, 
and this, reinforces the monocytic nature of microglia. 
For example Mill, a highly expressed HSC specific gene 
in the RMC, is a histone methyl transferase whose func- 
tional disruption is implicated in human leukemia [102]. 
Understanding the functions of these HSC specific genes 
may be important in comprehending the immune sys- 
tem related-roles of AMC and RMC. 

AMC express proliferation- and differentiation-related 
genes, Sox4, Soxll and Runxltl 

In order to validate the microarray data, we have 
analyzed the expression patterns of SOX genes (Sox4 
and Soxll) which are known to be involved in dif- 
ferentiation and Runxltl, which is involved in the 
proliferation of hematopoietic lineage cells. These 
genes were highly expressed by the AMC and their 
expression and role have not been studied in micro- 
glia, so far. Initially, nuclear expression of the tran- 
scription factor Soxll was shown to be associated 
with embryonic neurogenesis and lymphopoiesis 
[103,104]. However, there are no data on the role of 
SOX genes in microglia in which Soxll is expressed 
in the cytoplasm as reported in plasma myeloma 
cells and other B-cell lymphomas [104]. According 
to previous studies, the overexpression of the fusion 
protein, Runxl -Runxltl causes the downregulation 
of Csf-1 (a hematopoietic cytokine known to cause 
activation of microglia) [105] and Runx3 (a tumour 
suppressor) [106]. Our expression profile showed the 



increased expression of Runxltl in the AMC and 
downregulation of Csf-1 and Runx3 in AMC com- 
pared to RMC. Functional analysis of these tran- 
scription factors may help in understanding 
microglial proliferation and differentiation. 

Conclusions 

Overall, the transcriptome profiling has identified several 
genes, which help in elucidating morphological trans- 
formation and functions of AMC and RMC. These genes 
not only represent the physiological role of microglia in 
the developing brain but may also be useful therapeutic 
targets in neuropathologies in which microglia are 
implicated. 
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